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We studied numerically the Andreev scattering cross-sections of three-dimensional isolated quan- 
tized vortex rings in superfluid ^He-B at ultra-low temperatures. We calculated the dependence of 
the cross-section on the ring's size and on the angle between the beam of incident thermal quasi- 
particle excitations and the direction of the ring's motion. We also introduced, and investigated 
numerically, the cross-section averaged over all possible orientations of the vortex ring; such a cross- 
section may be particularly relevant for the analysis of experimental data. We also analyzed the role 
of screening effects for Andreev reflection of quasiparticles by systems of vortex rings. Using the 
results obtained for isolated rings we found that the screening factor for a system of unlinked rings 
depends strongly on the average radius of the vortex ring, and that the screening effects increase 
with decreasing the rings' size. 

PACS numbers: 

67. 30. he Vortices in He3 

67. 30. em Excitations in He3 

67.30.hb Hydrodynamics in He3 
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I. INTRODUCTION 

Supcrfluid turbulence consists of a disordered tangle of quantized vortex filaments which move under the velocity 
field of each other If the temperature, T is sufficiently smaller than the critical temperature, Tj,, then the normal 

fluid can be neglected and the vortices do not experience any friction effects Q • The simplicity of the vortex structures 
(discrete vortex lines) and the absence of dissipation mechanisms, such as mutual friction and viscosity, make superfluid 
turbulence a remarkable fluid system, particularly when compared to turbulence in ordinary fluids. Unfortunately, in 
contrast with ordinary turbulence, only few experimental techniques of flow visualization and detection of turbulent 
structures are available in superfluids at very low temperatures. Superfluid turbulence experiments are currently 
performed in both ^He and in ^He-B [oi-ll^. The methods used in these two liquids are different. In supcrfluid 
■^He-B, at temperatures T <C Tc, a powerful experimental technique, based on the Andrecv scattering of thermal 
quasiparticle excitations, can be used to detect the vortex filaments, see for example the review article [l5|. This 
technique, having been pioneered and developed at Lancaster University , is now also used at Aalto University in 

Helsinki [IJ] for measurements of vortex configurations. The Andreev scattering technique makes use of the fact that 
the energy dispersion curve, E = £'(p) of quasiparticle thermal excitations of momentum p is tied to the reference 
frame of the superfluid. From the Galilean invariance it follows that in this reference frame the dispersion curve 
tilts, becoming E{jp) + p • [l^, where Vg is the superfluid velocity. Thus, for thermal excitations whose energies 
are greater than the Fermi energy, ep (such excitations are known as quasiparticles), one side of the vortex filament 
presents a potential barrier and they are reflected back almost exactly, becoming quasiholes (excitations whose energy 
is smaller than ep)] the other side of the vortex lets the quasiparticles go through. Quasiholes are reflected or 
transmitted in the opposite way. The vortex thus casts a symmetric "Andreev" shadow for the quasiparticles at one 
side and for the quasiholes at other side, and by measuring the flux of excitations one detects the presence of the 
vortex. 

In our earlier works we developed a theory of ballistic propagation of thermal excitations near a single, rectilinear 
vortex filament in "^He-B [l6j . and studied interactions of thermal quasiparticles with simple, two-dimensional vortex 
configurations, such as clusters of vortex points and a gas of point vortices and/or vortex-antivortex pairs [13] . In 
the latter two works we found and investigated the phenomenon of the so-called 'partial screening' when the Andreev 
shadow of a system of vortices is no longer equal to the sum of shadows of individual vortices. However, the results 
following from our two-dimensional models should be regarded as qualitative rather than quantitative; a quantitative 
comparison with experimental observations and Andreev scattering data requires a fully three-dimensional study of 
vortex systems and quasiparticles trajectories. 

This work is concerned with the Andrecv reflection by individual vortex rings in three-dimensional geometry. Our 
study is particularly motivated by experimental observations of the transition from a gas of vortex rings to a dense 
vortex tangle and measurements of the decay of quantum turbulence generated by a vibrating grid shedding 
quantized vortex rings in alternating directions . Most conveniently the Andreev scattering of thermal excitation 
by quantized vortex rings can be characterized by the cross-section defined either by the ratio of the total number 
of quasiparticles reflected by the ring per unit time to the number flux density of quasiparticles incident on the ring, 
or, alternatively, by the cross-section defined as the ratio of the total power reflected by the ring to the flux density 
of energy carried by incident quasiparticles. In this work both cross-sections are calculated as functions of the ring's 
size and orientation with respect to the direction of the incoming beam of thermal excitations. Also calculated are 
the cross-sections averaged over all possible orientations of the ring. 

The plan of the paper is the following. In Sec. [Hi we shall introduce the equations of motion for ballistic quasiparticles 
in the superflow field, formulate the equations governing the fluid flow and the motion of quantized vortex rings, and 
define the cross-sections of interactions between thermal quasiparticles and vortex rings. In Sec. IIIII we describe the 
numerical method. In Sec. IIVI we shall calculate the scattering cross sections of the vortex rings and their systems. 
In Sec.|V]we shall draw the conclusions. 



II. BALLISTIC QUASIPARTICLES AND CROSS-SECTIONS OF THE ANDREEV SCATTERING IN 

THE FLOW FIELD OF QUANTIZED VORTEX RING 

We will be concerned with the propagation of thermal excitations in '^He-B at temperatures T <^ Tc, where 
Tc ~ I mK is the critical temperature. Below all numerical data are taken at bar pressure. 

Neglecting spatial variations of the order parameter, the energy of a thermal excitation of momentum p in the flow 
field Vs(r, t) generated by the quantized vortex ring is 



E{p, r, t) ^ y'ep2 + A2 + p.v,(r, <), 



(1) 
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where 

'v-7^.-^F (2) 

is the "kinetic" energy of a thermal excitation relative to the Fermi energy ep w 2.27 x 10~^^erg, p = |p|, m* w 
3.01 X m = 1.51 X 10^^'^ g is the effective mass of excitation in "^He-B (with m being the bare mass of the '^He 
atom) . We will be considering the propagation of thermal excitations at distances from the vortex core exceeding the 
zero-temperature coherence length, ~ 0.75 x 10^^ cm so that the superfluid energy gap can be regarded as constant, 
Ao = UQksTc « 2.43 x lO^^^erg (here fc^ is the Boltzmann's constant). Excitations with e^, > and < are 
called, respectively, quasiparticles and quasiholcs. 

Below we will follow the approach developed in our earlier works and assume that the interaction term, p-Vc, 

varies on a spatial scale which is larger than = ^Wi?/7rAo, where vp = \/'2.ep/m* « 5.48 x lO'^cm/s is the Fermi 
velocity. Then, following Refs. [T^H^, Eq. can be regarded as a semi-classical Hamiltonian for the excitation 
considered as a compact object (quasiparticle) , whose position and momentum are r(t) and p{t) respectively, yielding 
the equations of motion 



dE 



dp y^T+A^ 



(3) 



P = -^ = -TpP-v, , 4 
or or 

where a dot denotes a derivative with respect to time. Note that the right-hand-side of Eq. ([3|) represents the group 
velocity of thermal quasiparticle. 

In Eqs. ©-(HI), Vs represents the flow field generated by the quantized vortex ring. In the zero-temperature limit, 
the ring of radius R moves in the direction orthogonal to the ring's plane with the self-induced velocity (here we 
assume that the vortex core is hollow) 



2ttR 



In 



8R\ 1 
~^ ^ 2 



(5) 



where k = nh/m ~ 0.662 x 10~'^cm^/s is the quantum of circulation in '^He-B, and Qc is the core radius. Since Oc, 
being of the order of coherence length, is much smaller than the radius of the ring, it is appropriate to describe 
vortex lines as space curves of infinitesimal thickness. 

The details of the fluid velocity field, Vs(r, t) generated by the vortex ring self-propagating in the inviscid fluid can 
be readily found for example in monograph by Lamb plj . However, for the purpose of this study it will be more 
convenient, using periodic boundary conditions, to calculate the flow field numerically from the Biot-Savart law 

v.(r,0 = -|^/^xds, (6) 

where the integration extends over the whole vortex conflguration. The motion and evolution of a single vortex ring 
or a system of quantized vortices is governed by the equation 

|=v.(s,.), (7) 

where s = s{t) is a position of a point on the vortex line. 

We consider the Andreev scattering of the net flux of excitations which results in the case where there is a (small) 
temperature gradient. Assuming that the source of thermal excitations is sufficiently far from quantized vortices, 
the beam of quasiparticles incident on the vortex ring (or the vortex tangle) can be regarded as one-dimensional. 
The differential fiuxes of incident excitations, (nvdcg) (cm~^s~^), and energy, (nvdcgE) (ergcm~^s^^) (that is, 
respectiveh;, the number of quasiparticles passing and the total energy carried by these quasiparticles through unit 
area) are 

(nvdcg)^!^ N{E)vdcg{E)^^STdE, (8) 
(nvdcgE) [ N{E)vdcg{E)E^^6TdE , (9) 
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where 5T ^ T is a temperature difference between the source of excitations and the opposite side of the experimental 
cell, 

Np being the density of states at the Fermi energy with the corresponding Fermi momentum, 

Vdc, = ^ ^ > VF (11) 

is the group velocity of Bogoliubov quasiparticle, and f{E) is the Fermi distribution. At considered ultra-low tem- 
peratures, T < O.lSTc! typical of turbulence experiments in '^Hc-B, the Fermi distribution reduces to the Boltzmann 
distribution 

The quasiparticle trajectories resulting from interactions with the flow field of the vortex ring are determined from 
the solution of the problem represented by the closed system of equations ([3]), (O, and (I7|) (the details of the 
numerical method will be discussed below in Sec. The initial conditions follow from the Boltzmann distribution 
p2p and the assumption that the initial positions of incident quasiparticles are distributed randomly on the plane 
orthogonal to the beam of excitations. The solution of this problem yields the total number of quasiparticles Andreev- 
reflected by the vortex configuration per unit time, Nn (s~^), and the total power dissipated by Andreev- reflected 
quasiparticles, Qr (ergs~^). Then, the cross-section of Andreev scattering by the vortex configuration (ring) can be 
defined as either 

o-jv = NR/{nvdcg) (13) 

or 

(TE = QR/{nvdcgE) , (14) 

which we will call the particle and the thermal cross-section, respectively. Note that these cross-sections correspond 
to the area of Andreev shadow. Numerical calculations reported below in Sec. IIVI show that aN and as are practically 
indistinguishable in all considered situations, thus confirming that definitions (jl3p and (jl4p are correct. 

For the quantized vortex ring, the cross-section, a (below in this Section the subscript, iV or i? is omitted) is a 
function of the ring's size, i?, the ring's velocity, which is itself a function of R, and the angle a between the beam of 
incident quasiparticles and the direction of translational motion of the ring, i.e. cr = cr(i?, a). Different orientations 
of the ring with respect to the beam of excitations, and Andreev shadows in the cases where the ring moves either 
parallel or antiparallel to the direction of monochromatic beam of quasiparticles are illustrated on Fig. [1] 

In experiments (e.g. Refs. (lol. [ll|) vortex rings move in all possible directions. Of particular interest for interpre- 
tation of these and similar experiments will be the cross-section averaged over all possible angles a. Assuming equal 
probability for all ring's orientations, the probability that the ring's velocity is at an angle between a and a + da 
with the direction of the beam can be easily calculated as ^ sin a da, with < a < tt. Therefore, the angle-average 
cross-section of the ring of radius R should be calculated as 



1 r 

— / a{R, a) sin a da . 

2 Jo 



(15) 



III. NUMERICAL METHOD 



The superfluid velocity field, Vs(r, t) is calculated from Eqs. ©-([T]) by means of the vortex filament method using 
periodic boundary conditions. Calculations were performed in cubic periodic boxes of two sizes: a ~ 1.52 x 10~^ 
and a = 2.25 x 10~^ cm. The time evolution of the vortex configuration and, consequently, the velocity field are 
calculated by the second-order Adams-Bashforth method using the fixed time step At^ = 2 x 10~^ s. The technique 
of discretization of the vortex lines and the regularization of the Biot-Savart integral are standard. The details of our 
numerical algorithm are given in Ref. j2^ . 

Propagation of thermal excitations in the velocity field Vsi^, t) is governed by Eqs ©-([I]), which are solved using the 
numerical code based on the variable-step, variable-order implementation of the Numerical Differentiation Formulas 



(NDFs). A comprehensive and detailed description of the NDFs method is given in Ref. [23 1 
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FIG. 1. (Color online) a) Orientations of the quantized vortex ring with respect to the x-direction (shown by three red arrows) 
of the beam of quasiparticles. Andreev shadow (red dashed area) of the ring (blue solid line) moving b) parallel (in the positive 
x-direction, a — 0), and c) antiparallel (in the negative x-direction, a = tt) to the direction of the monochromatic beam of 
thermal excitations. 




FIG. 2. (Color online) The rectilinear line vortex (grey vertical line) and trajectories of the quasiparticles with average energy 
(E) = Ao + fcsT. The box shown in this figure is for visualization only - it is not the computational box; its sizes are 
-5 X 10"^ < a; < 5 X 10"^, -4 x 10"^ < t/ < 4 x 10"^, and -6 x 10"^ < 2 < 6 x 10"^ (units in cm). The purple (AB), 
blue (CD), green (EF), and navy (GH) lines are the excitation's trajectories corresponding to the initial impact parameters 
po — 0.0029, —0.0011, — 0.002 and — 0.003 cm respectively. The grey ellipse indicates the direction of the velocity field of 
the vortex, and the perpendicular to the vortex axis red (QF) line shows the extension {~ 0.002 cm) of the shadow casted by 
the vortex for the quasiparticles with energy (E). The solid and dashed arrows indicate the directions of motion of incident 
quasiparticles and retroreflected quasiholes, respectively. 

We tested our numerical method for Andreev reflection of quasiparticles, whose initial momentum is in the x- 
direction, by a single rectilinear vortex line located at a; = y = aligned along the z-direction; the velocity field of 
such a vortex is time independent. The energy of a thermal excitation, defined by Eq. ([T]), and its z-component of the 
orbital angular momentum, = PxV — PyX are both integrals of motion. Fig. [5] ihustrates trajectories of excitations 
with different starting conditions identified by the impact parameter po = yo, where yo is the initial y-coordinate 
of thermal excitation. (In other words, the impact parameter, po is defined as a minimum distance between the 
rectilinear trajectory, which the ballistic quasiparticle would have followed in the absence of a vortex, and the vortex 
core, see Fig. 2 in Ref. [iBl-) Fig [3] shows that in this calculation the relative errors in the quasiparticle's energy, 
AE/Eq = {E — Eq)/Eq and momentum, AJ/ J^o = {Jz — Jzo)/ Jzo (where Eq and J^q are, respectively, the initial 
energy and z-component of momentum), are less than 2.5 x 10~^ and 10"'^, respectively, so that our method conserves 
the integrals of motion very well. 

To test our numerical method in the three-dimensional case for moving rings or their systems, we note first that 
in such a case the energy and the z-projection of the orbital angular momentum are no longer integrals of motion. 
However, they nearly are because the time scales of the quasiparticle motion and of the fluid motion differ by at least 
an order of magnitude so that the flow flcld can be regarded as frozen, see the next paragraph. Our calculations show 
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FIG. 3. (Color online) a) Relative tolerance of the energy, {E — Eq)/Ei.). b) Relative tolerance of the ^-component of orbital 
angular momentum, [Jz — Joz)/Joz- Time is in units of /is. The red, black, blue and green lines correspond to the trajectories 
AB, GHG, CDC and EF of Fig. ^ respectively. 



that, in the three-dimensional case, for a single ring as well as for a system of several rings the relative errors in energy 
and in the orbital angular momentum remain sufficiently small in all considered situations, AE/ Eq < 1.5 x 10~^ and 
AJz/JzQ < 1.5 X 10^^, respectively. These estimates have been obtained for the smallest of considered rings, with 
R = Rmm = 7.5 X 10~^ cm. Clearly, the errors becomes smaller with increasing the ring's size. 

In principle, Eqs. ©-([Tl) and ©-Q should be regarded as a system of equations whose solution yields simultane- 
ously the time-dependent vortex configuration, the fluid velocity field, and the trajectories of quasiparticle thermal 
excitations. A numerical solution of such a system of equations presents formidable difficulties, in particular because 
the equations ©-(U]) governing the motion of quasiparticles are stiff. However, because a typical time of travel of an 
excitation within the computational box is much smaller than the characteristic timescale of the motion and evolution 
of quantized vortices, the problem can be reduced to simpler calculations of quasiparticle motion in the frozen flow 
field of vortex conflguration at any instant of time. To justify this approach we calculated, in the computational box 
of the larger size, a = 2.25 x 10"^ cm, the time of travel of the excitation in the velocity field of the rectilinear vortex 
illustrated in Fig. [2] We found that the longest average time, which is of the order of 10"'^ s, is spent within the box 
by quasiparticles whose initial momentum is p ~ pp -I- 2 x 10~^pp, where pp = sJIm^tF « 8.28 x lO^^gcm/s is 
the Fermi momentum. The average group velocity of these, slowest quasiparticles is about lO^cm/s. On the other 
hand, the largest velocity of the vortex points can be estimated as vg ^ k,/ (2'i:Rmin), where Rmin is the radius of the 
smallest vortex ring. In our calculations Rmin = 7.5 x 10~^ cm, so that vg ~ lOcm/s, which is an order of magnitude 
smaller than the average group velocity of the slowest quasiparticles. This justifies the 'frozen flow field' approach for 
solving Eqs ©-Q- 

In our numerical simulations the flux of thermal excitations is modelled by Nqg ~ 52272 quasiparticles entering 
from one side of the computational box and moving parallel to the x-dircction. Initial positions on the (y, 2:)-planc 
and energies of quasiparticles, Eq (Aq < Eq < 1.7Ao) are uniformly distributed. The quasiparticles arc characterized 
by the set of three integer numbers, (n, m, k), where n and m refer to the initial position of quasiparticle on the 
{y, z)-plane as follows: 



+ mSz„ 



n, m 



1, 



(16) 



a - a 
--+ndyn, Zm = -- 

whcrc N ~ 66, a is the size of the cube, 6yn and (5z„ are the distances, in y and z directions, respectively, between 
the nearest quasiparticles. The third number, k refers to the energy level corresponding to the discrete momentum, 
Pk = Pf + kSp, (where A; = 1, Nk) and is calculated as 



Enmk = \J {(-p)l + Al + Tpk ■ Vs(-a/2, y„, z„) 



(17) 



where {tp)k = Pfc/(2m*) — ep. In our calculations Nk ~ 12 (so that N x N x Nk = Nqs)- 

The incident number and energy differential fluxes of quasiparticles are now calculated numerically as follows: 



(nvdcg) 
(nvdcgE) 



n,ra,k 



{vdCg)nmk 

An 



mk 



,-17.6£„„fc/Ao 



(18) 
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FIG. 4. (Color online) The cross-sections (in units of cm^) of Andreev scattering by an isolated vortex ring, irjv (solid lines) 
and (j_E (dashed lines) as functions of the angle, a between the incident beam of excitations and the direction of motion of the 
vortex ring. The pairs of curves, from top to bottom, correspond to the ring's radii R = 3.9 x 10~^, 2.4 x 10~^, 1.4 x 10"'^, 
and 1.2 x 10^"* cm. 



where, using the definition of the gap parameter, Aq, for the considered temperature, T = O.lTc the exponent in the 
Boltzmann's factor has been replaced by — 17.6i?/Ao, and the numerical approximation for the quasiparticlc's group 
velocity, Vg — dE/dp follows from Eqs. ^ and ([T7| in the form 

t ^ \ ('^p)*: Pk , f a \ 

(vdCg)n,nk = —=====— + Vs --, ■ (19) 



VM^ + A^m* 

Having solved equations of motion ([3])-(|4]), the total number of quasiparticles Andreev-reflected per unit time, Nr, 
and the total power dissipated by Andreev-reflected quasiparticles, Qr, are calculated as 

Nr\ _ n {vdCg)niyik 



qD^T. Xo""" ( El, ) e-^'-'^""'=/^" A5„,„ , (20) 

n,7n,k 

where Rnmk = 1 if the (n, m, k)-th particle is Andreev-reflected, otherwise Rnmk = 0. In Eq. ((20|). ASnm = 
ji^y-n + Syn+i){Szr,i + Szm+i) IS the area element of the (y, z)-plane. 

As was mentioned earlier in this Section, the calculations of cross-sections have been performed in computational 
boxes of two different sizes: a = 1.52 x lO^^cm and a ~ 2.25 x 10~^cm. Reported in the next Section HVl the results 
of calculations of scattering cross-sections turn to be independent of the size of computational box. This justifies both 
the correctness of definitions (jl3p - (|14p and the accuracy of numerical approximations described in this Section. 



IV. RESULTS AND DISCUSSIONS 



We start with the calculation of the scattering cross section as a function of the radius of the quantized vortex ring, 
R, and the angle, a between the beam of incident quasiparticles and the direction of translational motion of the ring. 
The angle is a = in the case where quasiparticles and the ring move in the same direction, and a = tt if they move 
in opposite directions. 

The results of the calculation are shown in Fig. |4l As it has already been mentioned in Sec. [TTl the cross-sections 
ctat and ctb, defined by formulae (jl3|) and (|14p . respectively, practically coincide. This is hardly surprising considering 
that both of them correspond to the area where quasiparticles arc Andreev reflected, i.e. to the area of Andreev 
shadow. Unless specified otherwise, in the remainder of this Section we will not distinguish between and as, 
hence omitting the subscripts 'W and "£"'. From the results shown in Fig. |4]it is seen that the cross-sections (and 
hence the Andreev reflection area) is the largest in the case where the ring moves exactly towards the source of 
excitations (a = tt); in the case where the direction of the beam and that of the ring's motion coincide (a = 0), the 
cross-section is slightly smaller. The minimum reflection area occurs for angles slightly smaller than a = tt/2. For 
small rings the cross-section is almost angle-independent; most hkely, this is because at small intervortex distances 
the process of Andreev reflection becomes dominated by the partial screening effects investigated in our earlier works 

[Hill. 

Fig. [5] shows the dependence of the angle-average cross-section, defined in Sec. Ullby Eq- (US]), on the radius of the 



8 




R 

FIG. 5. (Color online) Angle- averaged particle (upper blue line) and thermal (lower green line) cross-sections (cm^) as functions 
of the radius of the ring (cm) . Inset shows the behavior of a with R for small rings. 



10 



<a> 

0.5 




0.05 0.1 

FIG. 6. (Color online) Scattering cross section (cm^) of a monochromatic beam of quasiparticles with energy E vs dimensionless 
parameter Se = (£ — Ao)/Ao for radii of the ring, from top to bottom, R = 3.34 x 10"^, R = 2.41 x lO"'*, and R = 4.83 x 10"'* cm. 

ring, R. For sufficiently large rings, R > 2.42 x 10~^cm the angle-average cross-section, (cr) exhibits almost linear 
dependence on R, which can be approximated as 

(a) w KR + C, (21) 

with K = 6.4 X 10^"^ cm and C = —5.13 x 10^^ cm^. For smaller rings, the behavior of (a) with R is also practically 
linear, with K = 4.3 x 10^'^ cm and C = 0. A rather sharp change of behavior occurring at i? « 2.42 x 10^^ cm can 
again be attributed to significant contribution of partial screening effects for small rings [IJl . 

So far we analyzed the Andreev scattering of quasiparticles whose initial energies are uniformly distributed in the 
incident beam. Of a certain interest are also cross-sections of a monochromatic beam, i.e. such that all incident 
quasiparticles have the same fixed energy, E. Illustrated by Fig. [SJ our numerical calculation shows a strong decrease 
of the angle-average cross-section with the non-dimensional parameter Se = (E — Ao)/Ao and the ring radius, R. 

Until now we have considered the Andreev reflection of quasiparticles on a single quantized vortex ring. Here we 
will analyze briefly the case where the incident beam of quasiparticles is Andreev-rcflccted by a system of n uninked 
quantized vortex rings. Such a system can also be characterized by the average radius of the ring, R = "Y^Ri/n. 
Enforcing the same total line length, Ltot for all n, we calculate the cross-section of Andreev scattering by the 
system of vortex rings for n progressively increasing from rii corresponding to a system of just a few large rings, 
to n = Uf corresponding to a system of many smaller rings, see Fig. [T] In our calculation rii = 6, Uf = 36, and 
Ltot = 0.908 x 10~^ cm. We have to emphasize that such a sequence of configurations of vortex rings is not due to 
the Biot-Savart evolution of the system, but is the result of numerically enforced algorithm. Our calculation shows 
that, due to the screening effects, the total cross-section of n rings, (7n{R) is smaller than the sum of angle-averaged 
cross-sections of individual rings, that is 

n 

a,,{R) < 5](a(i?,)) . (22) 

i=l 
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FIG. 7. (Color online) Modeling sequence of configurations of n unlinked vorex rings with the total line length being preserved 
for all n. Initial configuration of six rings (left), and final configurations of 36 rings (right). 
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FIG. 8. (Color online) Sum of the angle-average cross-sections of n individual vortex rings, n{a{R)} (upper green line, left), 
total (screened) cross-section, (t„ (lower blue line, left), and the screening factor (right) as functions of radius (cm). For each 
n all rings have the same radius, R = Ltot/{2Trn), with Ltot being the same for all n. Units of cross-sections are in cm^. 



Numerical calculations also show that the screening factor, defined as 



Sarel = 1 - (T„(i?) ^(tT(i?,)> , (23) 



increases with decreasing the average radius, R. For each n, fro = 6 to Uf = 36, we analyzed a number of 
configurations of the system of vortex rings, and found that for each n the total cross-section of the system oscillates 
around the value corresponding to the case where all rings have the same radius, R = R. This case is investigated in 
some more detail, and the calculated values of the total cross-section, cr„ , the sum of angle-averaged cross-sections of 
individual vortex rings, J^i'^i^i)) = ^lalR))^ and the screening factor, Sard = 1 — <^n{R) I {'n{a{R))) are represented 
in Fig. [5] as functions of radius, R = Lto(/(27rn). Results shown in Fig. |S] (left) indicate that, despite the total vortex 
line length remains the same for all values of i?, the total scattering cross-section of the system decreases substantially 
with radius. Figure [8] (right) shows the dramatic increase of the screening from 41% up to 69% with decreasing the 
rings radii from R ~ 2.41 x 10^^ cm to R = 0.402 x 10"'^ cm. This may be explained with the help of results illustrated 
by Fig. ini which indicate that the main contribution to the cross-sections of smaller rings is made by the low energy 
quasiparticles. For the high energy quasiparticles sufficiently small rings are almost transparent. When the number 
of rings is increased so that the rings' sizes are reduced, most of the low energy quasiparticles are reflected by the 
front-line rings, and just a small fraction of excitations reaches the rings in the bulk of the system; hence, because 
most of the high energy quasiparticles are not Andreev reflected at all, the screening effect increases. 
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V. CONCLUSIONS 

In conclusion, we have analyzed, for the first time, the three-dimensional Andreev reflection of thermal quasiparticle 
excitations by quantized vortex rings in "^He-B. The particle and thermal cross-sections (i.e. the Andreev reflection 
areas) of quantized vortex rings are defined and calculated; the results show a strong dependence of the cross-section 
on the angle between the incident beam of quasiparticles and the direction of motion of the vortex ring. It is also 
shown that the particle and the thermal cross-sections practically coincide. Of a primary interest for interpretation 
of experimental data is the cross-section averaged over all possible orientations of the vortex ring. This is calculated 
and its dependence on the size of vortex ring is analyzed in detail. It is apparent that the phenomenon of partial 
screening investigated in the authors' earlier works in two dimensions, plays a major role for rings of sufficiently small 
size in three dimensions. The results are generalized for the case of Andreev reflection by the system of vortex rings. 
It is found that due to the screening effects the total cross-section of the system of vortex rings is significantly smaller 
than the sum of cross-sections of individual vortices. Furthermore, were two system of vortex rings have the same 
total line length, the Andreev scattering cross-section is significantly larger of a system consisting of bigger rings. 
We introduced a screening factor of a system of vortex rings and showed that it decreases strongly with the average 
ring's radius. Our results may be helpful for inferring quantitative properties of ballistic vortex rings produced by 
a vibrating grid at its low velocities, as in the experiment reported by Bradley et al. (lo| . Our results can also be 
used for detecting the transition, observed in the cited experiment, from a gas of vortex rings to the dense vortex 
tangle (based on the two-dimensional model of vortex points, a qualitative analysis of change of the Andreev reflection 
coefficient during such a transition was given in our earlier paper [isj ; our new results may allow a more quantitative 
analysis). 
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